Mersenne Binomials and the Coefficients of the 
non- associative Exponential 

by L. Gerritzen 

Abstract 

The non-associative exponential series exp(x) is a power series with 
monomials from the magma M of finite, planar rooted trees. The coeffi- 
cient a(t) of exp(x) relative to a tree t of degree n is a rational number 
and it is shown that 



a(t) :-- 



a(t) 



2«- i -nr=i 1 (2 i -i) 



is an integer which is a product of Mersenne binomials. One obtains 
summation formulas 

^2&(t) = u(n) 

where the sum is extended over all trees t in M of degree n and 

nn-l 

n! ±A - 

The prime factorization of u(n) is described. The sequence (co(n)) n >i 
seems to be of interest. 



Introduction 

For a natural number n, we denote by M n the n-th Mersenne number 2 n — 1. 

n 

The Mersenne factorial n\u is denned to be the product Y\ Mi of all Mersenne 

i=i 

number M 1 , M n while the Mersenne binomial is defined to be 

r! A f(n - r)\ M 
for natural numbers r between and n. 
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In quantum calculus, see for instance [KC], there are the notions [n] for the q - 
analoque of n G N, [n]\ for the q - factorial and ["] for the q - binomial coefficients 
which are polynomials in the variable q. If /(2) denotes the value of a polynomial 
/ after substituting 2 for q, then [n](2) = M n , [n\\{2) = n\ M and [™] (2) = (™) M . 

Let M be the magma with unit 1m freely generated by a single element x. It can 
be identified with the set of all finite planar binary rooted trees, see [DG]. 

Let Q{{x}} the Q- algebra of power series with monomials from M. It was 
proved in [DG] that there is a unique series exp(x) G Q{{x}} such that exp(x) = 
1 + x+higher terms and 

exp(x) ■ exp{x) = exp{2x) 

see also Proposition (3.2). Let a(t) be the coefficient of the non-associative ex- 
ponential exp relative to the planar binary rooted tree t and put 

a(t) = n\ • uj{n) • a(t) 

where 

, v 2 n -\n-l)\ M 
W (n) := ^ . 

A main result of this article states that all a(t) are integers which are obtained 
as products of Mersenne binomials. The projection of the non-associative expo- 
nential to the classical one leads to the decomposition 

J^a(t) = u(n) 

if the sum is extended over all trees of degree n. 

In section 1 the notions of Mersenne order and Wieferich exponents at odd primes 
is defined. One gets results about the prime factorization of Mersenne numbers. 

In section 2 we show that all cu(n) are integers by computing the order ord p u(n) 
for all primes p. 

The applications about the coefficients a(t) of exp{x) are deduced in section 3. 
I would like to thank Doron Zeilberger for helpful hints. 



1 Mersenne orders 

We call M n = 2 n — 1 the n-th Mersenne number for n G N, n > 1. 

This is in contrast to the definition by D. Shanks, see [S], Chap. 1, which is also 

widely used in which M n is called Mersenne number only if n is a prime. 

Let p be an odd prime. 
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Definition 1.1. The Mersenne order v{p) at p is the smallest number r E N, r > 
1, such that p divides M r . 

As 2 V ^ = 1 mod p and t> (p) is the order of the class of 2 in (Z/pZ)*, we get that 
v(p) divides the order of the multiplicative group (Z/pZ)* which is (p — 1). 
Also p divides M n , if and only if n is a multiple of v(p). 
Also p < 2 V ^ from which follows that log 2 (p) < v(p). 
It has the consequence that lim p ^ 00 v{p) = oo. 

Definition 1.2. Let e{p) := ord p (M v ^) for any odd prime; it is called the 
Wieferich exponent at p. 

Always e(p) > 1 and p is called a Wieferich prime if e(p) > 2. 

Up to now only two Wieferich primes are known, namely 1093 and 3511 and it 

was checked that these are the only Wieferich primes < 4.10 12 in [CDP]. 

The Mersenne order at 1093 is 364 = ^y 2 - = 2 2 • 7 ■ 13 and the Mersenne order at 

3511 is 1755 = = 3 3 • 5 -13. 

The Wieferich exponent for both primes is equal to 2. 
Proposition 1.3. Let p be an odd prime, n E N, n > 1. 

(i) ord p (M n ) = 0, if n is not a multiple ofv(p). 

(ii) Ifn = v{p) -p r ■ k and p does not divide k, the ord p {n) = r and ord p (M n ) = 
e(p) + r 

Proof. Statement (i) is well-known. 
Let now n = v (p) ■ p r ■ k as in (ii). 

Let m be an integer and consider the group (Z/p m Z)* of units of the ring Z/p m Z. 

Let U rn be the subgroup of (Z/p m Z)* of class congruent to 1 modulo p. 

Then U m has order p m ~ l because any element a in U rn has a unique representation 

m— 1 

by a number a = 1 + ^^ with a,i E {0, ...,p—l}. The order of a is m—min({i : 
i=i 

aj 7^ 0}, m) because if a = 1 + p £ ■ x, x any number not divisible by p, then 

a p = l+p £+1 -y 

where 

y = % + EG>¥ £ ^ 1 

which shows that y = x mod p as p > 2. 

Let now a = 2^ where 2 is the class of 2 in Z/p m Z. 
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Then a G U m and a is represented by a number a = 1 + where x is not 

divisible by p. Then a pr is represented by a number 

1 + . y 

with y not divisible by p. 

It follows that ordpM n = e(p) + r because the order of a pr is equal to the order 
of a prk . □ 

Corollary 1.4. 

n 

where the product is extended over all odd primes p for which the Mersenne order 
v{p) divides n. 

2 Factorial Mersenne quotients 

F(n) := nr=i % is called the Mersenne factorial of n e N; it is also denoted by 
n\ M . 

Let 

, , 2"- 1 • F(n - 1) 

w « : = ; 

n\ 

for n > 1; it is called the factorial Mersenne quotient at n. Let m G N, m > 1, 
and p be a prime. There is a unique p-adic expansion 



m 

i=0 



^rriip 

with rrii G {0, 1, — 1}. Let 

r 

dp(m) := 



rrii 

8=0 



Proposition 2.1. ^ ord 2 Ld(n) = d 2 (n) — 1 
(nj If p is an odd prime, then 

(n — d p (n)) — (m — d p {m)) 



ord p (u(n)) — s(p) ■ m 



(P-1) 



where v(p) is the Mersenne order atp, m is the largest integer < (n—l)/v(p) 
and e{p) is the Wieferich exponent at p. 
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Proof. 1) m = d p (m) mod (p — 1) because m = ^^tth mod (p — 1). It is well 

i=0 

known that 

ord p (n\) = H - ^ = f 

A ] (P-1) ttV 

1) ord 2 (uj(n)) — (n — 1) — ord 2 (n^ because F(n) is odd. Thus 

ord2{u){n)) = d2(n) — 1 > 

for n > 1. 

2) Let now p be an odd prime. Then 

ord p F(n — 1) = ord p Mi. 
i=i 

If i = k ■ p r v{p) and A; is not divisible by p, then 

ordpMi — e{p) + r) 
and ordpMi = if p is not divisible by p. 

ord p (F(n - 1)) = ordpM jv{p) 
3=1 

where m p (n) = m := [^j] is the largest integer < (n — l)/t> (p). It follows 
from Proposition 1.3 that 

ord p M ja(p) = e{p) + ord p (j) 

as ordp(j) = ord p (jv(p)). As 



^ord p (j) = or-dp{m\) = — — ^ p | m ) 
we get the formulas 



t r^/ -, \ /x m ~ ^p(^) 
ora v b \n — \) — e(p) ■ m H — 

and 

(n - rfp(n)) - (m - d p (m)) 
ordp{uj{n) = e{p)m - 



p — 1 

□ 
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Proposition 2.2. uj(n) G N. 

Proof. We will show now that ord p uj{n) > 0. From the formula in 3) we get 

, . . .. (n — dJn)) — (m — dJm)) 

ordJuin)) >m- ± p --^ ^ 

p — 1 

> m _ n ~ d p^ > \ n ~ 1 i _ n z > o 

- [ P-1 J p-1 - 

because m — rf p (m) > 0, m > ^fy, n — d p (n) < n — 1 and "^j") ^ 

(3) From 1) and 2) we get that ord p (uj(n — 1)) > for all primes p. As oo(n — 1) 
is a rational number, it follows that u(n — 1) is an integer. 

□ 

Example 2.3. (1) ord 7 uj(13) = 3, because n = 13, m = ^y- = 4, d 7 (4) = 
4, d 7 (13) = 7, e(7) = 1, n-d 7 (n) = 6, m-d 7 {m) = ord 7 (cu(13)) = 4-| = 3 

(2) ord 7 cj(100) = 21, because v(7) = 3,n = 100, m = [^] = 33, d 7 = 
100) = 4,d 7 (33) = 9 ord 7 {uj{\m)) = 33 - 5^ = 21, because log(u(n)) « 
( n + 1 )%(2) + n-(n + l/2)/o( 7 (n). 

Let 7Tm(^) be defined to be number of odd primes p such that the Mersenne order 
v(p) is < x — 1. 

It seems interesting to determine the asymptotic behaviour of ir M . 

One can expect some results about ttm from the prime factorization of ti\m, 

because 7Tm(^) is the number of primes dividing n\u- 

Example 2.4. 7Tm(16) = 15 because the primes dividing 16!m are the odd primes 
< 16 and the primes 17, 23, 31, 43, 73, 89, 127, 151, 257, 8191. 

3 Coefficients of the non-associative exponen- 
tial 

Let M be the magma with neutral element 1m freely generated by a set consisting 
of one element x. The multiplication in M is a map • : M x M — > M and the 
restriction of this map onto (M — {1m}) x (Mx{1 m }) is injective and its image 
is M — {1m, x}. It means that for any t G M — {1m, x} there is a unique pair 
(ti, t 2 ) G M x M,ti ^ 1m t 1 h, such that t = h ■ t 2 . 

There is a unique homomorphism deg : M — > N such that deg{lM) = 0, deg(x) = 
1. 

It is important to realize that any t G M, t ^ 1 M , gives rise to a unique planar 
binary rooted tree with deg(t) leaves, see [1], section 1, p. 163. The grafting of 
such trees corresponds to the multiplication map of M. 
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Let P = Q{{x}} be the Q-algebra of power series with monomials from M. Thus 
any / G P has a unique expansion 

/ = £<#)•* 

with c(t) G Q. It will also be called the algebra of tree power series in x or the 
algebra of power series over Q in a nonassociative and noncommutative variable 
x. 

Define the order ord(f) of / to be min{deg(t) : c(t) 7^ 0}, if / 7^ and ord(0) = 
+00. The following proposition is obtained by standard calculus methods. 

Proposition 3.1. (i) There is a unique Q-linear map j- : P — > P such that: 



and 



^ - 1 



for all f,geP. 

f := 4-(f) is called the derivative of f with respect to x. 

(ii) Let g G P with ord(g) > 1. Then there is a unique Q-algebra homomor- 
phism rj g : P — > P such that r) g (x) = g. One denotes r) g (x) also by f o g or 
f(g) and calls i] g the substitution homomorphism induced by g. 

(Hi) There is a canonical Q-algebra homomorphism [ ] : P — > Q[[x]], where 
Q[[x]] denotes the classical Q-algebra of power series in an associative (and 
commutative) variable x such that [x] = x. 

One calls [/] the classical power series associated to f G P. 
Proposition 3.2. There is a unique f G Q{{rc}} such that 

f (0) = 1 
/(*) • /(*) = f(2x) 

Moreover f'(x) = f(x) where f'(x) denotes the derivative of f relative to x. The 
tree power series f(x) is called the nonassociative, noncommutative exponential 
series and will be denoted in this article by e x or exp(x). 

Proof. Inductively we define a map a : M — > Q by putting a(t) :— 1, if t — 1m 
or t = x and otherwise 

ogi) • a(t 2 ) 
w 2 n - 2 
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if n = deg(t) and t = t\ • t 2 with U e M. 

Put / := Y^a(t)-t E P = Q{{x}}. It is easy to check that f(x) ■ f(x) = f(2x) : 
teM 

f(2x) = J2 a (t) ■ 2deg(t) ■ 1 
teM 

f{x) ■ f{x) = ^i)ait 2 )h ■ t 2 

(ti,t2)eMxM 

= 1 + 2a(x) + 

+ (a(t) ■ a(x°) + a(x°)a(t) + a(t!)a(t 2 )) ■ t. 

teM,deg(t)>2,t=(ti,t 2 ) 

As 2a(t) + a(ti) ■ a(t 2 ) = 2 n a(t) for any f 6 M with deg(t) > 2,t = h ■ t 2 , we 
obtain the desired functional equation. 

For the proof that f'(x) = f(x) we refer to [DG]. □ 

From the proof we get that the coefficients a(t) of exp satisfy: a(l) = a(x) = 1 
and 

q(ti)a(t 2 ) 

where n = deg{ti ■ t 2 ). 
Let a(t) := n • a;(n) • a(t). 

Proposition 3.3. Let t±,t 2 G M,ni = deg(ti) > 1 and t = t\ ■ t 2 ,n = deg(t) = 
rii + n 2 . Then 

n-2 

ait) = i — —)Ma(t 1 )a(t 2 ) 



Proof. 



i / \ a(h)a{t 2 ) 
ait) = n\ ■ u{n) 



= m 



2 n - 2 
u(n) a(ti) d(t 2 ) 



2M n _i n 2 !cj(n 2 ) 
2 n " 1 F(n-l) a(ti)a(t 2 ) 

2M n _i 2™i- 1 F(m - l)2 n 2" 1 F(n 2 - 1) 
Fin -2) 



Fin, - l)F(n 2 - 1) 
n-2 



a(ti)a(t 2 ) 



n± — 1 



a(ti)a(t 2 ). 

M 



□ 
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It follows from (3.3) that a(t) is a product of Mersenne binomials. More precisely: 
Corollary 3.4. a(t) = TT I n ^ a \ ^ ] e N w/iere /(t) the set of inner 

nodes of t and n(a) is the degree of t< a where t< a is the tree below a which is 
defined to consist of all nodes b for which the simple path from b to the root of 
t is passing through a. Also n\{a) is the degree of the left factor s\ of t< a (such 
that t< a = si • s 2 ). 

Example 3.5. Let T n be the set of trees in M of degree n defined as follows: 
T\ = {x}. 

If T n _i is already defined, then T n = x ■ T n _i U T n _i • x. 

Then %T n = 2 n ~ 2 for n > 2. One can show that a(t) = 1 iff t e T n . 

Corollary 3.6. 

a(t) = iv(n) 

deg(t)=n 

Proof. We consider the canonical algebra homomorphism 

{] : ®{{x}} ^ Ql[x]]. 

OO „ 
^ ^£ 

Then [exp] is the classical series ~- As [t] = x n if t <E M has degree n, it 

n=o n ' 

follows that a(t) — — . □ 



71! 
deg(tj=n 



For 1 < k < n - 1 let 



deg(t-± ) — k 
deg(t2)=n—k 



Then 



Sk{n)=( U k _ 2 \ «(*i)«(*2) 

\ / M deg(ti)=k 

deg(t'2) — rL — k 

Sk(n) = ^ _ 2 ^j ■ u(k) ■ uj(n - k) 

n-l 

It follows from Corollary 3.6 that u(n) = J2Sk{n). Thus 



k=i 

n-l 



fn - 2\ 

Corollary 3.7. > u(k) ■ u(n — k) = iv(n). 

k=1 \ k - l /M 



Example 3.8. u(l) = uj{2) = 1, u(3) = 2, w(4) = 7, w(5) = 42 = 2 • 3 • 7 
Thus 

W (6) = w (*M 6 " fc ) 

fe=l \ / M 

= 1 • w(5) + M 4 • 1 • w(4) + M \' M3 co(3) 2 + M 4 • w(4) • 1 + 1 • u(5) 
= 2-42 + 30 -7 + 5 -7-4 = 434 = 2- 7-31 

How to understand from this partition of u(6) that 31 is a factor of u>(6)? 
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